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ABSTRACT 

We present a determination of the optical QSO luminosity function and its cos- 
mological evolution with redshift for a sample of over 6000 QSOs identified primarily 
from the first observations of the 2dF QSO Redshift Survey (2QZ). For QSOs with 
—26 < Mb < —23 and 0.35 < z < 2.3, we find that pure luminosity evolution (PLE) 
models provide an acceptable fit to the observed redshift dependence of the luminosity 
function. The luminosity function is best fit by a two-power-law function of the form 
$(Lb) oc [(Lb/Lb) 01 + {Lb/L'qYY 1 . Exponential luminosity evolution models, both 
as a function of look-back time, L B (z) — L B (0)e klT , and as a general second-order 
polynomial, L B (z) oc \Q k i z + k 2z ^ were f oun d to provide acceptable fits to the dataset 
comprising the 2QZ and the Large Bright Quasar Survey. Exponential evolution with 
look-back time is prefered for qo — 0.05, while the polynomial evolution model is 
prefered for q = 0.5. The shape and evolution of the LF at low redshifts (z < 0.5) 
and/or high luminosities, not currently well sampled by the 2dF QSO survey, may 
show departures from pure luminosity evolution, but the results presented here show 
that over a significant range of redshift, PLE is a good description of QSO evolution. 

Key words: galaxies: active - quasars: general 



1 INTRODUCTION 

The QSO optical luminosity function (OLF) and its evolu- 
tion with redshift provides fundamental information on the 
overall demographics of the AGN population. It provides 
constraints on the physical models for QSOs (Haehnelt & 
Rees 1993; Terlevich & Boyle 1993), information on models 
of structure formation in the early Universe (Efstathiou & 
Rees 1988) and a picture of the ionizing UV/optical lumi- 
nosity density from QSOs as a function of redshift (Meiksin 
& Madau 1993; Boyle & Terlevich 1998). 

Based on the major ultra-violet excess (UVX) QSO sur- 
veys of the 1980s (Bracessi: Marshall et al. 1983, Palomar- 
Green: Green, Schmidt & Liebert 1986, Durham-AAT: 
Boyle et al. 1990), a picture emerged in which the low- 
intermediate redshift [z < 2.2) QSO OLF was modelled by 
a two-power-law function with a steep bright end (4>(Lb) oc 
L- 3 ' 6±01 ) and a much flatter faint end ($(L B ) oc L~ U2±0 - 1 ) 
whose redshift dependence was best fit by pure luminosity 
evolution, i.e., a uniform increase in luminosity toward high 
redshift (see e.g., Marshall 1985; Boyle et al. 1988; Hartwick 



& Schade 1990). This evolution was modelled as a power- 
law in redshift of the functional form L B oc (1 + z) klj , where 
3 < < 3.4 for an Einstein-de Sitter universe. 

Later QSO surveys, particularly those that focussed on 
bright magnitudes (LBQS: Hewett, Foltz & Chaffee 1993, 
EQS: Goldschmidt & Miller 1998, HBQS: La Franca & Cris- 
tiani 1997) have reported evidence for a more complex form 
of evolution in which the bright end of the OLF showed a 
significant steepening with increasing redshift. Indeed the 
steepening was sufficiently dramatic that, at the very lowest 
redshifts studied, the OLF showed a single featureless power 
law (Koehler et al. 1997) with little evidence for any cosmic 
evolution amongst the brightest QSOs (Mb < —27). Cou- 
pled with this, further claims were made (Hawkins & Veron 
1995) that the observed break in the luminosity function at 
higher redshifts was much less dramatic than had previously 
been reported. 

With the increasing number of high redshift (z > 2) 
QSOs discovered in surveys with well-defined selection cri- 
teria, evidence was also found that the strong power-law 
evolution does not continue on beyond z ~ 2. The nature 
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of this change was modelled in a variety of ways; from a 
simple halt to constant comoving density (see Boyle et al. 
1991) between 2 ~ 2 to z ~ 3 to a slackening off of the 
evolution rate over the broad redshift range 1.6 < z < 3.3 
(Hewett et al. 1993; Warren, Hewett & Osmer 1994). From 
surveys of higher redshift QSOs (3.5 < z < 5.0), strong ev- 
idence emerged for the onset of a dramatic decline in the 
QSO space density at z > 3.5 (see e.g Warren et al. 1994). 

However, compared to our knowledge of the galaxy lu- 
minosity function, QSO OLF determinations were based on 
relatively few objects. This was particularly true for the ex- 
tremes of the redshift distribution z < 0.3 and 2 > 4 or the 
high luminosity end of the OLF (Mb < —27), where depar- 
tures from the pure luminosity evolution model were first 
noted. In these regimes, the numbers of objects in suitable 
surveys could be counted in the tens. Even in the most well- 
sampled regions of the QSO (Mb, z) plane, surveys com- 
prised a few hundred or up to one thousand QSOs in total. 

The long-heralded arrival of the new generation of large 
QSO surveys compiled with facilities such as the 2dF (Lewis 
et al. 1998) or Sloan Digital Sky Survey (Loveday et al. 1998) 
will shortly generate many tens of thousands of QSOs with 
which to carry out the definitive study of the QSO OLF and 
its evolution with redshift. These surveys also provide the 
opportunity to study the QSO OLF as a function of various 
physical parameters (e.g., emission-line strength, continuum 
slope) or types (e.g., BAL QSOs). 

In this paper, we present a new determination of the 
QSO OLF based on the first ~ 6000 QSOs identified in the 
2dF QSO redshift survey (2QZ) (sec Smith et al. 1998). The 
2QZ is currently more than a factor of 10 larger than pre- 
vious QSO surveys to a similar magnitude limit (bj < 20.85 
mag). When complete, it is planned that the 2QZ will com- 
prise 25000-30000 QSOs. In section 2 we present a brief 
overview of the survey and in section 3 we describe the anal- 
ysis of the survey. We present our conclusions in section 4. 



2 DATA 

2.1 The 2dF QSO redshift survey 

2.1.1 Input catalogue 

For the purposes of the analysis we have used the current 
version (as of September 1999) of the 2QZ catalogue con- 
taining 6684 QSOs in total. The identification of the QSO 
candidates for the 2QZ was based on broadband ubjr colour 
selection from APM measurements of UK Schmidt (UKST) 
photographic plates. The survey area comprises 30 UKST 
fields, arranged in two 75° x 5° declination strips centred 
on 8 — —30° and S = 0°. The 8 = —30° strip extends from 
a — 21 h 40 to a = 3 h 15 and the equatorial strip from a 
= 9 h 50 to a = 14^50. The total survey area is 740 deg 2 , 
when allowance is made for regions of sky excised around 
bright stars. The 2QZ area forms an exact subset of the 
2dF Galaxy Redshift Survey (GRS: see Colless 1998) area, 
with identical 'holes' used for both surveys. In a typical 2dF 
field, approximately 225 fibres are devoted to galaxies, 125 
to QSOs and 25-30 fibres are devoted to sky. The data are 
reduced using the standard 2dF pipeline reduction system 
(Bailey & Glazebrook 1999). 

In each UKST field, APM measurements of one bj plate, 
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Figure 1. Coverage function for the 2dF QSO redshift survey. 
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Figure 2. Photometric completeness contours in the (6j, 2) plane 
for the 2QZ. The contours are plotted at 5 per cent intervals: 
dotted lines denote 5, 15, 25, ...95 per cent completeness; solid 
lines correspond to 10, 20, 30. ..90 per cent. 



one r plate and up to four u plates/films were used to gener- 
ate a catalogue of stellar objects with fej < 20.85. A sophis- 
ticated procedure was devised to ensure catalogue homo- 
geneity (Smith 1998). Corrections were made for vignetting 
and field effects due to variable desensitization in the UKST 
plates, these effects being particularly noticeable at the 
edges of plates. The criteria for inclusion in the catalogue 
were (it - bj) < 0.36; (u - bj) < 0.12 - 0.8(6j - r); (6j - r) < 
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0.05. Based on the colours of QSOs which have previously 
been identified in the survey region, we estimate the cat- 
alogue is ~ 90% complete for z < 2 QSOs and comprises 
~55% QSOs (see Croom 1997). Subsequent spectroscopic 
observations with 2dF have confirmed this QSO fraction, 
with the principal contamination arising from galactic sub- 
dwarfs and compact blue galaxies. Full details of the con- 
struction of the input catalogue may be found in Croom 
(1997) and Smith (1998). 

2.1.2 Spectroscopic observations 

2dF spectroscopic observations for the 2QZ began in Jan- 
uary 1997, although the bulk of the redshifts have been 
obtained in the observing runs after October 1998 as the 
2dF system has gained functionality (increased number of 
fibres, faster field re-configuration times). Over 6000 QSO 
redshifts have been obtained with the 2dF, and the 2QZ 
is now the largest single homogeneous QSO catalogue in 
existence. Each 2dF field in the survey is observed for typ- 
ically one hour. The spectral resolution is 4A pixel -1 and 
the spectra cover the wavelength range 3700A-8000A. This 
set-up gives a typical signal-to-noise ratio of approximately 
seven or greater between 4000A-6500A in the continua of 
the faintest objects (bj = 20.85) in the 2QZ. This was suffi- 
cient to identify up to 85 per cent of the objects in each 2dF 
field down to this magnitude limit. Where poor conditions 
prevented us from achieving this identification rate in any 
single 2dF observation, we included in the final catalogue 
only those QSOs that were brighter than the magnitude at 
which 80 per cent of the QSO candidates in the field had a 
positive identification. 



2.1.3 Spectroscopic Incompleteness 

In the OLF analysis below, we correct for incompleteness in 
the 2QZ by using an effective area which is a function of 
both apparent fej magnitude and redshift. Hereinafter, all fej 
magnitudes are corrected for galactic extinction using the 
Schlegel, Finkbeiner & Davis (1998) values. 

The 2QZ is subject to three forms of spectroscopic in- 
completeness. First, a small fraction of survey fields were ob- 
served in poor conditions, and do not reach the target spec- 
troscopic completeness of 80 per cent. As described above, 
the magnitude limits for these fields were made brighter un- 
til the desired completeness level was reached. Secondly, we 
corrected the inevitable trend to increasing spectroscopic in- 
completeness at fainter magnitudes by correcting the actual 
area observed by the fraction of objects with a reliable spec- 
troscopic identification. This was done in 0.025mag bins to 
track accurately the magnitude-dependent incompleteness. 

Finally, we applied a uniform correction factor of 0.71 
to the effective area to account for the 29 per cent objects 
in the input catalogue that have not yet had been observed 
by 2dF, despite being located in areas already observed in 
the 2dF survey. This is because the numbers of candidates 
in the combined QSO/galaxy catalogues in a typical 2dF 
field is ~ 50 per cent higher than the number of fibres avail- 
able. Near-complete (> 95 per cent) spectroscopic coverage 
of the survey area by 2dF is achieved by a complex tiling 
algorithm, with each 2dF field being visited 1.5 times on 



average. Since the observations are still at a relatively early 
stage, the coverage of the survey area with 2dF is relatively 
patchy and the completeness fraction is still relatively low 
(71 per cent). The resulting effective area of the survey is 
plotted in Figure [j]. 

2.1.4 Photometric Incompleteness 

Photometric incompleteness, arising from the errors in the 
photographic magnitudes (±0.1 mag in each band) and vari- 
ability (due to the noncontemparaneous nature of the ubjr 
plates on each field) will cause QSOs to exhibit u — bj/bj — r 
colours outside our selection criteria. The incompleteness is 
a complex function of both magnitude and redshift. An esti- 
mate of this incompleteness has been made by Croom et al. 
(in preparation) using mean colour-redshift relations (zero- 
pointed to the 2QZ system) from the non-colour-selected 
QSO survey of Hawkins & Veron (1995). At 2 < 2.3 the 
Hawkins & Veron colours accurately trace the mean colours 
of the 2QZ QSOs. The dispersion in the colours as a function 
of magnitude was derived from the observed dispersion of 
colours for the 2QZ QSOs in the redshift interval 1 < 2 < 2, 
the regime of highest completeness. A Monte Carlo sim- 
ulation using the measured dispersion in colours and the 
Hawkins & Veron mean colours, with 10 6 QSOs in each 
A&j and Az bin, was then used to predict the photometric 
completeness as a function of both magnitude and redshift. 
The completeness contours are shown in Figure ^. The pho- 
tometric completeness is largely independent of magnitude 
and is at least 85 per cent or greater over the redshift range 
0.4 < 2 < 2.1. At higher redshifts the completeness rapidly 
drops, falling to below 50 per cent at z > 2.3. We have there- 
fore chosen this redshift as the upper limit to our analysis 
below. Although the catalogue contains many hundreds of 
QSOs with z > 2.3 (which will be used for the clustering 
analysis that is not so dependent on photometric complete- 
ness) , any small errors in the completeness estimates at these 
redshifts can lead to large variations in the computed effec- 
tive area. 

At the lowest redshifts, the derived photometric com- 
pleteness is still relatively high (~ 60 per cent at z = 0.2). 
However, low redshift, low luminosity QSOs (Mb > —23) 
dominated by their host galaxy light will be lost from the 
2QZ as a result of the stellar selection criterion applied to the 
input catalogue. Although it may be possible to correct for 
such selection effects, for the present analysis we have sim- 
ply chosen to exclude all low luminosity QSOs (Mb > —23) 
from the OLF fitting procedure below. Given the bright ap- 
parent magnitude limit of the 2QZ (bj = 18.25 mag), this 
effectively limits the minimum redshift in the complete 2QZ 
sample to 2 > 0.35. We have therefore adopted this as our 
low redshift limit for all QSO surveys used in the OLF analy- 
sis below. Thus, out of a total of 6684 QSOs identified in the 
2QZ (13552 objects observed in 219 2dF fields), 5067 fulfil 
the criteria for inclusion in the complete sample (qo = 0.5, 
H = 50kms- 1 Mpc- 1 ). 

2.2 Other Samples 

We are currently in the process of extending the bright limit 
of the 2QZ from bj — 18.25 mag to bj = 17 mag using ob- 
servations made with the FLAIR spectrograph on the UK 
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Figure 3. (Mb, z) distribution for QSOs used in this analysis. Mb calculated using go = 0.5, Ho = 50km s 1 Mpc . The symbols used 
for the different surveys are indicated on the plot. The bright limit for the 2QZ survey, bj > 18.25 mag, is denoted by the solid line. 



Table 1. Parameters for the surveys used in our analysis. Nqso is the number of QSOs with Mb < — 23 (qo = 0.5, Hq = 50km s 1 Mpc x ) 
in the redshift interval 0.35 < z < 2.3. 



Survey 


Nqso 


Bum 


Area 


Reference 










(deg 2 ) 




2dF 


5057 


18.25 < 6j 


< 20.85 


196 


This paper 


LBQS 


867 


16.50 < bj 


< 18.85 


454 


Hewett et al. (1995) 


EQS 


106 


15.00 < 6j 


< 18.00 


330 


Miller et al. (unpublished) 


HEQS 


23 


15.00 < 6j 


< 17.65 


611 


Koehler et al. (1997) 


PG 


36 


13.00 < bj 


< 16.67 


10653 


Green et al. (1986) 



Schmidt Telescope. These observations are not yet complete 
and we have incorporated a number of independent QSO 
surveys at brighter magnitudes to extend the present analy- 
sis of the OLF to higher luminosities. Table [l] lists the areas, 
magnitude limits, and the number of QSOs within the com- 
pleteness limits (M B < -23, 0.35 < z < 2.3) for each QSO 
survey used in this paper. The absolute magnitudes and red- 
shifts for the QSOs in these surveys are plotted in Figure [| 

The Large Bright QSO Survey (LBQS: Hewett et al. 
1995) provides a complementary sample to the 2QZ. The 
details and completeness of the survey are well established 
and it provides a large number of QSOs in the 1.5 — 2 mag 
interval brighter than the bj > 18.25 2QZ survey limit. 

At the brightest magnitudes B < 16.5 mag, the largest 
currently published survey is the Palomar-Green survey 
(Green et al. 1986). However, the completeness of this survey 



has been called into question by a number of recent surveys 
including the Edinburgh Quasar Survey (EQS: Goldschmidt 
et al. 1998) and the Hamburg/ESO Quasar Survey (HEQS: 
Koehler et al. 1997). Unfortunately, full details of both these 
catalogues have yet to be published, and although we have 
access to the unpublished EQS, details of its completeness 
have yet to be accurately determined. 

Given the remaining uncertainty over their details, and 
apparent discrepancy between QSO surface densities derived 
from these brighter surveys, we chose to carry out our anal- 
ysis of the OLF for three separate survey combinations: a) 
2QZ + LBQS; b) 2QZ + LBQS + PG; c) 2QZ + LBQS + 
EQS + HEQS. Given the available information, the 2QZ + 
LBQS data sample constitutes our primary data sample in 
this paper. When more extensive published catalogues are 
available from brighter surveys (in particular the HEQS) 
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they will clearly provide a powerful test of the OLF models 
at the brightest absolute magnitudes (see Figure ||). 

A further survey at intermediate magnitudes, the Ho- 
mogeneous Bright QSO Survey (HBQS) is currently under 
construction (see Cristiani et al. 1995). There is consider- 
able overlap between the survey areas for the HBQS and 
EQS. The currently published HBQS samples contains fewer 
QSOs than the EQS sample to which we have access. We 
have therefore chosen not to include the HBQS in the anal- 
ysis below. 

With one exception, the surveys are all largely inde- 
pendent of one another. The the current coverage of the 
2QZ results in a very small overlap with the LBQS; there 
is only one QSO in common between the 2QZ and LBQS. 
Given the complexity of the current 2dF coverage function 
in the spatial domain (due to incomplete tiling) , we therefore 
chose to treat the LBQS and 2QZ as independant surveys. 
At brighter magnitudes, the PG, HEQS and EQS are all 
independant of one another. 

The only exception is the overlap between the LBQS 
and EQS. There are 32 QSOs in the EQS sample which are 
also in the LBQS sample, or approximately 30 per cent of the 
EQS sample used in this analysis. The QSOs are distrubuted 
throughout the four UK Schmidt fields in common between 
the two surveys. In the analysis below we have therefore 
excluded these four fields from the EQS during the model 
fitting procedure, but we test the acceptibilty of the best-fit 
model against the full EQS sample. The number of QSOs 
for the EQS quoted in Table 1 refers to the full sample. 



3 ANALYSIS 
3.1 1/V estimator 

As the first step in our analysis, we obtained a graphical 
representation of the OLF and its evolution with redshift. 
We used the new 1 /V estimator devised by Page & Carrera 
(1999) to derive a binned estimate of the OLF. Although it 
addresses one of the potential biases of the traditional 1/V a 
statistic (Schmidt 1968), like most binned LF estimators 
the Page & Carrera (1999) method assumes that absolute 
magnitude and redshift bins can be chosen sufficently small 
so that the effects of evolution and a steeply rising OLF 
are negligable across each bin. For previous samples this 
has been difficult to achieve whilst still retaining a sufficient 
number of QSOs in each bin. 

With the large numbers of QSOs now available for this 
analysis, we were able to minimise these effects by choosing 
much smaller bins in redshift and absolute magnitude. We 
used 0.25 mag bins and 10 redshift bins equally spaced in 
log(l + z) over the interval 0.35 < z < 2.3 to compute the 
1/V estimate of the OLF. 

For an OLF that evolves significantly across a Mb ,z bin, 
1/V methods will still yield biased estimates of the OLF in 
bins cut by the magnitude limit of the sample. In such bins, 
the estimated OLF will only represent the mean space den- 
sity of objects sampled in bright Mb , low z region of the bin. 
In the absence of any magnitude limit most objects would 
tend to lie in the faint Mb , high z part of the bin (at least for 
a power law LF undergoing strong redshift evolution) and so 
the estimate of the OLF in bins which are not fully sampled 



will be biased low. The extent of the bias will, of course, de- 
pend on the extent to which the OLF evolves across the bin. 
Our choice of small bins alleviates, but does not remove this 
bias and so we have also excluded all Mb , z bins which con- 
tain the bj ~ 20.8 mag survey limit in our binned estimate 
of the OLF. A similar bias also affects (but in the opposite 
sense) bins cut by the bright mangitude limit of the sam- 
ple. However, we also imposed the constraint that we would 
not compute the OLF in bins where there were five or fewer 
objects. This latter conditions prevents any bins cut by the 
bright survey limit from appearing in estimate of the OLF 
below. 

The resulting OLF, calculated for a flat universe with 
go = 0.5 and a Ho = 50km s _1 Mpc -1 , in plotted Figure ^. 
Absolute Mb magnitudes were derived for the QSOs using 
the k-corrections derived by Cristiani & Vio (1990). The 
median number of QSOs in each plotted bin is 40, although 
some bins contain up to 170 QSOs. The shape and redshift 
dependence of the OLF in Figure ^ are strongly suggestive 
of a luminosity evolution model similar to those derived pre- 
viously from the Durham/ A AT sample (Boyle et al. 1988). 



3.2 Maximum Likelihood Analysis 

3.2.1 Method 

To obtain a more quantitative descriptions of suitable mod- 
els, we carried out a maximum likelihood fitting procedure 
for a number of models to the data (see Boyle et al. 1988). 
This technique relies on minimizing the likelihood function S 
corresponding to the Poisson probability distribution func- 
tion for both model and data (Marshall et al. 1983). We 
tested the goodness-of-fit of the model to the data using 
the 2D Kolmogorov- Smirnoff (KS) statistic. The KS is no- 
toriously insensitive to discrepancies between the data and 
the model predictions in the wings of the distributions. This 
problem is particularly severe in the present analysis, where 
the weakness of the KS test at the brightest absolute mag- 
nitudes is compounded by the vanishingly small fraction of 
objects in this region of the Mb,z plane in the combined 
sample (see Figure ^). 

To alleviate this problem, we derived 2D KS statistics 
from each survey separately; combining them into a final 
KS probability using the Z statisitic described by Peacock 
(1983). This approach enables each of the model predic- 
tions to be tested against the data in specific regions of the 
(Mb, z) plane sampled by different surveys. A significant re- 
jection of the model, even by a relatively small sample, will 
thus have major impact on the overall acceptability of the 
fit. 

In the fitting procedure, we used no more free param- 
eters in any model than was required to obtain an accept- 
able fit. We defined a priori an acceptable fit as one which 
could not be rejected at the 99 per cent confidence level or 
greater i.e., a KS probability (Pks) of 1 per cent or greater. 
Errors on the fit parameters correspond to the AS* = 1 con- 
tours around each parameter, or, equivalently, the 68 per 
cent confidence contour for one interesting parameter. 
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Figure 4. Luminosity function for the 2QZ + LBQS dataset in a flat go = 0.5 universe 



3.2.2 Models 

Guided by the appearance of the OLF in Fig ^ we chose to 
model the OLF $(!/, z) as a two-power4aw in luminosity ,Q 

$ ( l b z) = 

1 ' ' [(L s /L* B )" + (L B /L* B n 
Expressed in magnitudes this becomes 



All simpler forms of the OLF e.g., single power law, Schechter 
function were strongly ruled out (Pks < 10 — 10 ) in preliminary 
fits. 



#(M B ,«) = 



10 0.4[(a+l)(M B -Af*(*))] + 10 0.4[(^+l)(M B -M*(z))] ' 



where the evolution is given by the redshift dependence of 
the break luminosity L B or magnitude, M B (z). 

We first attempted to fit the evolution by using 'stan- 
dard' models, i.e., an exponential luminosity evolution with 
fractional look-back time (r), such that: 

Lb{z) = L B (0)exp(kir), 

and a power-law luminosity evolution model with a redshift 
cut-off: 
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Table 2. Best-fit OLF model parameters from maximum likelihood analysis. 



Surveys sampled 


Evolution 
ivioaei 


go 


Nqso 


Q! 


P 


\ i * 
B 






Mpc ^maj 


r -i 






2dF & LBQS 


exp(kr) 


0.05 


6100 


3.37 


1.55 


-21.16 


7.11 


— 


0.47x10" 


6 


0.52x10" 


l 




expikr) 


0.50 


5924 


3.32 


1.39 


-19.67 


6.98 





0.15x10" 


5 


0.14x10" 


3 




A- 12 - k 2 z 2 


0.05 


6100 


3.29 


1.47 


-22.49 


1.30 


-0.25 


0.56x10" 


6 


0.24x10" 


2 




k\z — k 2 z 2 


0.50 


5924 


3.43 


1.58 


-21.99 


1.35 


-0.27 


0.11x10" 


5 


0.12 




2dF, LBQS & PG 


expikr) 


0.05 


6136 


3.37 


1.63 


-21.33 


6.98 




0.41x10" 


6 


0.14x10" 


1 




expikr) 


0.50 


5960 


3.52 


1.57 


-20.00 


6.95 




0.10x10" 


5 


0.32x10" 


6 




k±z + k 2 z 2 


0.05 


6136 


3.30 


1.55 


-22.33 


1.42 


-0.29 


0.52x10" 


6 


0.34x10" 


2 




k\z + k2Z 2 


0.50 


5960 


3.60 


1.77 


-22.39 


1.31 


-0.25 


0.68x10" 


6 


0.15x10" 


1 


2dF,LBQS, HEQS, k, EQS 


expikr) 


0.05 


6265 


3.95 


1.87 


-21.80 


7.16 




0.17x10" 
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In contrast to previous analyses, the power-law evolu- 
tion model provided a poor fit to all datasets (Pks <S 0.01) 
for all cosmological models. This implied that any decline in 
the evolution at high redshifts (z ~ 2) was not well repre- 
sented by an abrupt cut-off at z max - We therefore adopted 
a more general exponential evolution law incorporating a 
second-order polynomial function (hereinafter referred to as 
polynomial evolution) that gives a smoother transition in 
the power law behaviour of the evolution at z > 2 as allows 
for the possibility of negative evolution at high redshift Q 

k 1 z + k 2 z 2 



L B (z) =Lb(0)10' 
or equivalently, 

M B 0) = M B (0) - 2.5(fei2 + k 2 z 2 ). 
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3.2.3 Results 

The best-fit parameters, and KS probabilities for the polyno- 
mial and exponential evolution fits are given in Table ^. The 
statistical errors on individual parameters are Aa ~ ±0.05, 
A/3 ~ ±0.1, AM B ~ ±0.2, Afci (exponential) ~ ±0.1, 
Afci (polynomial) ~ ±0.05, Afc 2 ~ ±0.02. 

Acceptable fits (Pks > 0.01) to the primary 2QZ + 
LBQS sample were found for both evolution models. The ex- 
ponential model favoured low go and the polynomial model 
favoured high go. 

The predicted differential number-magnitude, n(m), 
and number-redshift, n(z), relations for the 2QZ survey 
based on the best-fit go = 0.5 polynomial evolution model 
to the 2QZ ± LBQS samples are shown in Figure |H[ The 
model predictions provide a good fit to the derived n(m) 
and observed n(z) relations for QSOs with 0.35 < z < 2.3 
and Mb < —23 from the 2QZ survey also plotted in this 
figure. 

The extrapolated maximum L B in the polynomial evo- 
lution model occurs in the range 2.46 < z < 2.50 for 

t Polynomial evolution models were first fit to the radio 
galaxy/QSO LFs by Dunlop k. Peacock (1990) 



' 

0.2 0.4 0.6 0.8 1 

Look— back Lirne t 

Figure 6. Evolution of L* B as a function of fractional look-back 
time (r) for the polynomial evolution model in a go = 0.5 uni- 
verse. The solid line denotes the best-fit evolution model obtained 
from fitting to the 2QZ + LBQS data-sets. The dashed line rep- 
resents the extrapolated evolution of L* B between z = 2.3 and 

2 = 5. 



0.05 < go < 0.5. The behaviour of L B as a function of frac- 
tional look-back time for the go = 0.5 universe is shown in 
Figure |(| 

In Figure |^ we have also plotted the Poisson signifi- 
cance of the residual difference between the best-fit go = 0.5 
polynomial evolution model and 2QZ + LBQS data-set for 
the Mb,z bins used in the 1/V analysis of the OLF. It 
can be seen that, within the range of luminosity-redshift pa- 
rameter space covered by the data, there are no significant 
(> 3a) differences between the model and data. Although 
there may be a possible weak trend at low redshifts for the 
model to over-predict the numbers of QSOs at low luminosi- 
ties and to under-predict numbers at high luminosities, any 
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Figure 5. Upper panel: Derived differential number-magnitude, n(m), relation for the 2QZ survey (filled dots) and prediction from 
best-fit polynomial evolution model (qo = 0.5) to LBQS+2QZ sample. Lower panel: Observed and predicted number-redshift, n(z), 
relation for the 2QZ survey. Model as for upper panel. Note that the derived n(m) relation from the 2QZ survey has been corrected for 
both spectroscopic and photometric completeness, whereas the n(z) relation corresponds simply to the observed numbers of QSOs in the 
2QZ survey. Correspondingly the predicted n(m) relation comes directly from the model, whereas the n(z) model prediction has been 
corrected for the assumed survey incompleteness. 



such discrepancies are not yet statistically significant with 
this current dataset. 

We found no significant difference in the values of the 
best-fit parameters when the fits were restricted to z > 0.5 
or 2 < 2.0, i.e., excluding the sample regions where we apply 
the largest correction for incompleteness. The extrapolated 
peak in occurs at the same redshift (z = 2.5) for model 
fits to both the z < 2 and z < 2.3 samples see (Figure []). 

The addition of the bright samples yielded much poorer 
overall fits. Inclusion of the PG sample gave, at best, 
marginally acceptable fits for the exponential evolution 
model (go = 0.05) and polynomial evolution model (go = 
0.5). No acceptable fit was found for a go = 0.5 universe 
with the incorportion of HEQS and EQS, and only a barely 
acceptable fit (Pks = 0.01) was obtained for the qo = 0.05 
exponential evolution model. 



Previous attempts to characterize departures from pure 
luminosity evolution (PLE) have led to models with a 
redshift-dependent bright-end slope (La Franca & Cristiani 
1997). We attempted to fit such a model for the polyno- 
mial evolution models to the data-set including the HEQS 
and EQS samples using a simple redshift dependence of the 
form 

a(z) = q(0) + n-iZ. 

For qo = 0.5 the best fit value of K3 = 0.36, while for 
go = 0.05, K3 = 0.0. However, this model did not provide 
an acceptable fit to the data-set (Pks < 0.01) for either 
go = 0.5 or go = 0.05. 

By virtue of its size, rejection of the fits by the brighter 
data-sets was dominated by the EQS. Nevertheless the EQS 
contains relatively few QSOs, and a full understanding of 
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Figure 7. Greyscale respresentation of the P oisson significance 
of the discrepancy ((N Cl \ 3s — AT pIe< j)/ y/ N pYe ^) between the ob- 
served number of QSOs (N ^, s ) in the 2QZ + LBQS data-set and 
predicted number of QSOs (N prc ^) from qo = 0.5 polynomial evo- 
lution model. The Mg , z bins chosen are the same as those used 
in the 1/V determination of the OLF. Only bins with A r prod > 5 
have been plotted. 

the evolution of the very brightest QSOs clearly still re- 
quires both larger samples and a more detailed knowledge 
of their properties. This should be available with the com- 
pleted HEQS. 

There could be a variety of reasons why luminous QSOs 
may exhibit deparatures from pure luminosity evolution 
models. For example, bright QSO samples are known to con- 
tain a greater fraction of radio-loud QSOs (Peacock, Miller 
& Longair 1986) and also contain a greater fraction of gravi- 
tationally lensed QSOs (Kochanek 1991). Both effects could 
give rise to systematic departures from simple luminosity 
evolution models fit predominantly to less luminous QSOs. 
Equally, these results may imply that luminous QSOs simply 
evolve differently from the bulk of 'normal' QSOs. 

We also explored the effect of different cosmological 
models on the OLF for our primary 2QZ + LBQS sam- 
ple. The results are reported in Table pi For an Q,a = 
universe we fit the polynomial evolution model to the data 
for a variety of different values of go- All values of go in the 
range < go < 0.5 yielded acceptable fits, with probabilty 
reaching a broad maximum in the the region go = 0.3 — 0.4. 
However, this value is highly dependent on the evolution 
model chosen to fit the data. For example the exponential 
evolution model favours lower values of go (go < 0.1, see 
Table 2). It is unlikely that the QSO OLF can be used to 
define useful contraints on the value of go until a meaningful 
physical model for QSO evolution is available. 

The best-fit OLF parameters for a non-zero cosmologi- 
cal constant in a flat universe (i.e., f2 m + = 1) were also 
obtained using the expressions for comoving distance r(z) 
and comoving volume element dV/dz given by: 



^(z) = r(z)V«m(l + z) 3 + n A 

A variety of fits for different values of Qm and Q,a are pre- 
sented in Table ^. Luminosity evolution remains a good fit 
to the 2QZ + LBQS dataset in all cases. As above, it is ap- 
parent the OLF provides little potential for discrimination 
between cosmological models with a non-zero Qa- 



4 CONCLUSIONS 

For absolute magnitudes —26 < Mb < —23 (go = 0.5) and 
redshifts 0.35 < z < 2.3, pure luminosity evolution (PLE) 
has been shown to produce an acceptable fit to the QSO 
distribution. This region of parameter space contains the 
QSOs reponsible for the vast majority of the AGN luminos- 
ity density in the redshift range between 0.35 < z < 2.3. The 
best-fitting PLE models exhibit an exponential increase in 
luminosity with either look-back time (go = 0.05) or as a 
second-order polynomial function of redshift (go = 0.5). In 
the near future, the space density of z > 4 QSOs predicted 
from extrapolations of these models may be usefully com- 
pared with the significant numbers of such QSOs now being 
identified at these redshifts by the Sloan Digital Sky Survey 
(Fan et al. 1999). 

Inclusion of bright surveys such as the EQS and HEQS 
do result in departures from luminosity evolution similar in 
form to those seen by other authors (e.g. Hewett et al. 1993, 
La Franca & Cristaini 1997, Goldschmidt & Miller 1998). 
The principal region of parameter space in which such de- 
partures are seen lies at z < 0.5 where those authors claim 
a significant flattening of the luminosity function. The 2QZ 
contains relatively few QSOs with 0.35 < z < 0.5 and does 
not yet probe the redshift range z < 0.35. Hence there is at 
present no inconsistency between the various studies. Our 
analysis of the 2dF survey does however show that, over the 
redshift range 0.35 < z < 2.3 and for absolute magnitude 
Mb < —23, pure luminosity evolution does provide an ac- 
curate phenomenological model of QSO evolution. 
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